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The generation of ELF (of the order of 10 Hz) currents and electromagnetic fields in the 
off-equatorial E-region (90-120 km) of the Earth's ionosphere and their subsequent pene- 
tration into the deeper ionospheric layers is studied theoretically and numerically. These 
ELF currents and fields are generated by the interaction of an electromagnetic pulse with 
the E-region at its lower boundary located at 90 km above the Earth's surface. The wave 
penetration (with a typical wavelength of the order of 10 km) of the generated ELF currents 
and fields into the deeper ionospheric layers up to 120 km takes place due to the dominance 
of the Hall conductivity over the Pederson conductivity in the region between 90-120 km and 
penetration becomes diffusive above 120 km. During night time, the increase in the wave 
speed due to the reduced conductivities leads to the deeper penetration. As the angle be- 
tween Earth's magnetic field and horizontal is increased (going away from the equator) , the 
currents and fields penetrate deeper into the ionospheric layers with increased wavelength, 
the magnitudes of horizontal (east-west) and vertical currents decrease near the boundary 
and the vertical electric field decreases drastically. The horizontal (east- west) current inte- 
grated along vertical is fitted with a current distribution which can be replaced by a line 
current raised above its actual height by the half width of the current distribution for the 
purpose of the calculation of its radiation. The maximum of the total east-west current (310 
Amps) remains same for various simulation parameters due to the magnetic shielding. 
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I. INTRODUCTION 

The short scale size (50-100 km) electromagnetic fields in the ELF range are generated below 
Earth's ionosphere by variety of phenomena such as lightning discharges [3-0], impulsive fields 
created by seismic events [7i] or by ground based horizontal electric dipole antennas [8]. These 
fields can interact with the E-region (90-120 km) of the Earth's ionosphere and generate currents 
and associated electromagnetic fields which will radiate in the presence of the conducting Earth. 
Infact, inspired by the physics of the quasi-static equatorial electrojet (EEJ) a novel concept of 
the conversion of an ELF ground based horizontal electric dipole antenna (HED) into a vertical 
electric dipole (VED) antenna in the E-region was proposed [8]. At ELF frequencies, the VED 
radiates 10 5 times more efficiently than HED of the same dipole current moment 0]. In the case 
of quasi-static EEJ, tidal motions drive and maintain horizontal (zonal) electric fields of the order 
.5-1 mV/m perpendicular to the ambient magnetic field over long times (several minutes to hours) 
and over a 600 km strip in the E-region of the dip equatorial ionosphere. This electric field drives 
downward Hall current. At steady state and to zero order, current continuity requires that a 
vertical polarization electric field be built to prevent the downward Hall current from flowing. This 
electric field is larger than the zonal electric field by the ratio of the Hall-to-Pedersen conductivity, 
approximately a factor 30, resulting in vertical electric fields in excess of 10 mV/m and associated 
predominantly with eastern electrojet currents of more than 10 A/km. As noted by Forbes [1981], 
current continuity requires the presence of a vertical current, with current closure established by 
field aligned currents. These currents result in ground based quasistationary magnetic fields of 100 
nT or more. 

The interaction of quasistatic electric fields with the equatorial E-region leading to the genera- 



tion of EEJ has been studied extensively both experimentally and theoretically However 



the interaction of ELF range electromagnetic fields with i?-region is relatively less understood. 



Eliasson and Papadopoulos 15f| studied generation and penetration of the ELF current and elec- 
tromagnetic fields into the equatorial E-region. They found that the interaction of the ELF pulsed 
and continuous wave fields with the equatorial E-region leads to the generation of both vertical 
and horizontal currents which penetrate into ionospheric layers as Helicon waves. It is the objec- 
tive of this paper to study the physics of the interaction of the ELF electromagnetic field with 
off-equatorial E-region where the Earth's magnetic field makes finite angle 6 with the horizontal. 
It is found that as the angle between Earth's magnetic field and horizontal is increased (going 
away from the equator), the currents and fields penetrate deeper into the ionospheric layers with 
increased wavelength, the magnitudes of horizontal (east-west) and vertical currents decrease near 
the boundary and the vertical electric field decreases drastically. 

The paper is organized as follows. The next section contains physics of the E-region, numerical 
model, boundary conditions and simulation setup. Numerical fits for the real ionospheric conduc- 
tivities used in the simulations are provided in this section. Section [TTT1 presents modelling results 
for various simulation parameters, viz., day-night time conditions, pulsed and continuous wave an- 
tenna field and various values of the angle between the Earth's magnetic field and the horizontal. 
In section HVl we summarize our results. 



II. INTERACTION MODEL AND SIMULATION SETUP 

The E-region of the Earth's ionosphere (90 to 120 km above the Earth surface) consists of 
partially ionized plasma magnetized by the Earth's magnetic field. In this region, ions are viscously 
coupled to neutrals through collisions (fj n >> w c j) while electrons are strongly magnetized (u en « 
(jj C e)- Here z^ n and v en represent ion- neutral and electron- neutral collision frequencies, and u c i and 
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oo ce represent ion and electron cyclotron frequencies respectively. The presence of the background 
Earth's magnetic field Bo drives currents perpendicular to the magnetic field, in addition to the 
parallel currents. This gives rise to the tensor conductivity of the ionosphere. The current is related 
to the electric field by generalized Ohm's law. 

E x B 

J = o-||E|| + o-pE_ L -o- H — — -, (1) 

-DO 

where an, op and an are parallel, Pederson and Hall conductivities respectively. The current 
parallel to the magnetic field Bo is controlled by an, along the perpendicular electric field Ej_ 
(Pederson current) by ap and perpendicular to both Ej_ and Bo (Hall current) by an- These 
conductivities depend on density, collision frequency and magnetic field, and are given by, 

°\\ = £ o | — + — (2) 

"en u m I 

•■in- I Venice v in^d \ ,., 

°P = e — ( 2 2 + 2 2 ) (3) 
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Here it has been assumed that the plasma consists of electrons and one ion plasma species. In 
the E-region, where ui ce » v en and uj c i « fj n , the dominant Hall conductivity gives rise to 
the Helicon waves which mainly govern the plasma dynamics. While in F-region (above 120 km), 
where uo ce » v en and uo c i » fj n , the Hall conductivity decreases and plasma dynamics is mainly 
diffusive. In our numerical modeling, we have approximated the real vertical profiles of day time 



conductivities ly] by numerical fits given by the following functions. 



1 (5) 

" exp(-z/L li ||) + a-2,1 1 exp(-z/L 2 ,||) 

1 (6) 

ai t p exp(-z/Li i p) + a 2 ,p exp(z/L 2i p) 

aH = 1 

ai tH exp(-z/ Li }H ) + a 2 ,p exp(z/ L 2 ,h) 

The values of the parameters a's and L's are listed in TableUand the profiles for day time conditions 
are plotted in Fig. [lb. For night time conditions we assume that the conductivities are decreased 
by a factor of five due to the decreased electron and ion densities. It can be seen that in E- 
region, Hall and Pederson conductivities increase with altitude and peak around 110 km and 130 
km respectively. The Hall conductivity dominates Pederson conductivity in the E-region while 
Pederson conductivity is dominant above 120 km. 



01,11 = 3.07 x 10 9 fhn 


Li ii = 5.36 km 


o 2 ,|| = 0.63 x 10 3 f7m 


L 2j || = 18 km 


a hP = 7.67 x 10 12 fhn 


Li p = 5.36 km 


a 2 ,p = 0.990m 


Li.p = 18.8 km 


a hH = 1.53 x 10 n fhn 


L\,h = 5.36 km 


a 2 ,H = 0.0157fim 


L 2 ,h = 10.67 km 



TABLE I. Parameter values used in the numerical fits ©-(JT]) of the day time conductivity profiles. 
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FIG. 1. (a) The geometry of the simulation. The ionospheric layer is above z — 90 km and free space 
is below z — 90 km. The simulation box extends horizontally from x — —3000 km to x — 3000 km and 
vertically from z = 90 km to z = 160 km. The constant geomagnetic held Bo makes an angle 9 with the 
positive x direction. The ELF magnetic field B a „t which couples to the ionosphere at z = 90 km is generated 
by an antenna placed at x = z = 0. (b) Numerical fits of crn, ap and <jh as given by Forbes and Lindzen 
[16j for day time conditions. For night time conditions, the conductivities are decreased by a factor of 5. 
(After 0.) 



The geometry of the simulation is shown in Fig. [T^,. The ionospheric layer is treated as a 
two-dimensional layer varying in the horizontal direction x and the vertical direction z, and no 
variation along y. In this model, free space is below z = zq = 90 km, while the ionospheric layer 
extends vertically above z = 90 km. The simulation box covers the region from z = 90 km to 
z = 160 km and from x = —3000 km to x = 3000 km. The ionospheric layer and free space are 
magnetized by a constant external geomagnetic field Bo which makes an angle 6 with the positive 
x direction. 

Similar to the previous work (la |. we assume a very simple model for the antenna, that of an 
equivalent infinite length wire in the y direction located at x = z = 0, so that, 

Rant(x, Z, t) = B ant (t) 9 _ ? ■ (8) 
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where B an t{t) is the value of the antenna magnetic field at the bottom of the E-region (z = zq) 
and x = 0, and the normalized coordinates are x = x/zQ,y = y/zo, and z = z/zq. 

We have chosen two kinds of time dependence for B ant (t), namely, pulsed antenna field and 
continuous wave antenna field. The time dependence of the antenna field for the two cases are shown 
in Fig. [2j For the pulsed case we use an antenna field of the form B ant = i?o,ant exp[— (t — to) 2 /2D 2 ] 
for t < t , B ant = Bo :ant for t < t < ii and B ant = B 0>ant exp[-(t - t 1 ) 2 /2D%] for t > h, 
where the maximum amplitude Bo^ an t = 1 nT is reached at time to = 0.05 s, and the pulse 
switched off smoothly at t\ = 0.15 s, using the pulse rise and decay time Dt = 0.01 s. For 
the continuous wave case we choose a 10 Hz antenna field that is ramped up smoothly so that 
Bant = B 0tant exp[-(t - tO) 2 /2D 2 ] sin(207rt) for t < t and B ant = B 0tant sin(207rt) for t > t , where 
the pulse reaches its maximum amplitude -Bo, ant = 1 n T at t$ = 0.5 s, and the rise time is Dt = 0.15 
s. 

It is convenient for numerical purpose to express the generalized Ohm's law in terms of the 
impedance tensor. 

E = pJ (9) 

where p is the impedance tensor which is obtained by inverting the conductivity tensor. For the 
background magnetic field making an angle 9 with positive x-direction, it is given by, 



p\\ cos 2 9 + pp sin 2 i 

— ph sin 9 
G°i| — Pp) sin 9 cos / 



Ph sin I 
PP 



(PI I 



- Pp ) sin f/ cos i 
Ph cos 9 



-pn cos 9 pn sin 2 9 + pp cos 2 i 



Pn Pl2 Pis 

—P\2 P22 P23 
Pl3 —P23 P33 



where /?|| = 1/cn, pp = crp/(ap + a 2 H ) and pn = c#/(<7p + & H ). The other equations governing 
the dynamics in the E-layer are Ampere's law and Farady's law. 



V x B 
dB 
~dt 

Combining Ampere's and Faraday's law with Ohm's law, we obtain an evolution equation for B. 



(i J 
-V 



E 



(10) 
(11) 



<9B I 

— = V x \p.V x B] 

dt po 



(12) 



Equation (|12p describes the Helicon wave dynamics when impedance tensor is nondiagonal and 
diffusive phenomena when impedance tensor is diagonal. Separating equation (|12p into components 
and using the condition V.B = 0, we get three coupled equations in B x , B y and B z . 
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Here V.B = is used in place of the x-component of the evolution equation (|12p to calculate B x . 
For 9 = 0°, equations (fl~3l) - ([T5l) and impedance tensor p reduce to the forms which were studied 
earlier by Eliasson and Papadopoulos [151 ] . 
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FIG. 2. The time dependence of the antenna magnetic field B ant (t) used in the simulations, (a) pulsed 
antenna field with pulse width 0.1 sec. (b) continuous wave antenna field with frequency 10 Hz. 



As boundary condition at the plasma- free space boundary z = zq, we use the continuity of the 
magnetic field B z and its z-derivative la ]. This also gives continuity of the parallel (to the boundary 
surface in the x-y plane) electric field. A detailed discussion about the boundary conditions is given 
by |l5]. 



III. RESULTS 

We have conducted a series of numerical studies of the system (|13p - (|15p with various parameters 
of interest, viz., day and night time conductivities, pulsed and continuous wave antenna field and 
different values of 9. In what follows we first present results for a fixed value of = 5° and 
vary other parameters. Then we change values of for all other parameters in order to see the 
^-dependence of the results. Finally, a distribution will be fitted to z-integrated J y profile, which 
is useful for the calculation of the radiation from J y . 

Figs. [3] and H] show spatial profiles of magnetic field, currents and electric fields for pulsed 
antenna field for the case 9 = 5° for day and night time conductivities, respectively. The left 
column shows results at t = 0.1 s when the antenna field is at its maximum (1 nT), and right 
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FIG. 3. The magnetic field components B Xl By and B z (nT), the current density components j x -, jy and j z 
(nA/m 2 ) and the electric field components E x , E y and E z (mV/m) for day time conditions, pulsed antenna 
field and 8 = 5° at t = 0.1 sec. (left column) and t = 1 sec. (right column) 



column in the relaxation phase at t = 1 s when the antenna field is zero. It can be seen at 
t = 0.1 s that interaction of the antenna field with the ionospheric plasma at the lower boundary 
at z = 90 km generates currents and associated electromagnetic fields structures. The oblique wave 
penetration with decreasing wavelength of these structures into the deeper ionospheric layers up 
to z ~ 110 km can also be seen. The x-scale ~ 100 km of the structures is larger than the z-scale 
~ 10 km. The amplitude of the vertical electric field E z ~ .6 mV/m is one order of magnitude 
larger than those of the horizontal components E x and E y . The vertical current J z has upward 
and downward flow structure and its amplitude is 1-2 order of magnitude smaller than those of 
the horizontal currents J x and J y . The disparity between the magnitudes of J x and J z , and x- 
and ^-scales is consistent with the current continuity condition dJ x /dx + dJ z /dz = 0. Another 
interesting feature to note here is that at x = and just below z = 90 km the horizontal magnetic 
field component B x ~ 2 nT which is almost double of the antenna magnetic field (1 nT), while 
vertical magnetic field component B z is much smaller than the antenna field. This happens due 
to the shielding of plasma interior from the incident magnetic field. At the lower boundary, the 
antenna electric field E y drives Pederson current J py = opE y along E y and Hall current Jh± m 
the x-z plane perpendicular to the background oblique magnetic field Bo = Bq(xcos9 + zsinO). 
A perpendicular electric field Ej_ develops to limit the magnitude of the Jh± to the value required 
by the current continuity condition. This drives Hall current Jh y = (JhE± in the negative 



8 



Might, R«_il-£ 



1 0O 





t s • 




□ 
□ 
□ 
□ 
□ 
□ 
□ 
□ 
□ 



□ OS 
0.04 

0.02 
O 

-D.D2 
0.05 
O 

-a_as 

0.02 



-0.02 
5 



Trf -1 13 

120 
10 



FIG. 4. The magnetic field components B Xl B y and B z (nT), the current density components j x , j y and 
j z (nA/m 2 ) and the electric field components E x , E y and E z (mV/m) for night time conditions, pulsed 
antenna field and 9 — 5° at t — 0.1 sec. (left column) and t = 1 sec. (right column) 



y-direction. The total current J y = Jh y + J py generates the magnetic field component B x which 
cancels the antenna field inside the plasma while at the ionospheric boundary at z = 90 km and 
x = 0, it adds to the antenna field to almost double its value. The z-component of the antenna 
magnetic field is partially cancelled by the magnetic field produced by the induced current J y . 

Due to the current continuity, the perpendicular current is shunted parallel to Bo- The z- 
component of this parallel current and of the Hall current J#j_ give rise to the upward-downward 
flow structure of J z . The current components in terms of electric field components can be written 
as, 

Jx = | cos 2 8 + ap sin 2 6)E X — an sin 8E y + (an — op) sin 9 cos 9E Z (16) 
Jy = &h sin 6 Ex + crpEy — ajj cos 9E Z (17) 
Jz = (^"11 — 0"p) sin 9 cos 9E X + an cos 9E y + (an sin 2 9 + ap cos 2 9)E Z (18) 

In the region 90 km < z < 110 km, a\\/a H ~ 50 - 150 >> 1 and a\\/a P ~ 2400 - 1900 >> 1. 
Therefore for small values of 9 (for example for 9 = 5°,cos# ~ 1 and sin# ~ 0.08) J x is mainly 
dominated by parallel current as E x ~ E y ~ O.l^, while J z has contributions both from parallel 
and Hall currents in x-z plane. For large values of 9 close to 90°, J x and J z exchange their roles. 
The out-of-plane current J y is mainly Hall dominated as an » ap for z < 120 km and E z >> E y . 
The phase speed of the wave penetration can be estimated from whistler dispersion relation 
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= oj/k = fcii/9ff//io- This gives = 100 km/s for an average value of pn = 4000 S/m in the 
region 90 km < z < 110 km and a half wavelength of 100 km giving k\\ = 7r/100 = 0.314 km -1 . 
This gives penetration distance ~ 10 km by t = 0.1 sec., as is seen in Fig. [3l From whistler 
dispersion relation, A oc pu ~ for given values of ui and fc||. Since an increases from z = 90 

km to z = 110 km, the wavelength decreases, which is consistent with the simulations. 

In the relaxation phase (t = 1 sec.) when antenna field has decayed to zero, the sign of the 
structures have reversed and the magnitudes have decreased. The sign is reversed because in 
the decay phase of the incident pulse the antenna electric field changes its sign, which generates 
currents and electromagnetic fields in the opposite direction. The waves penetrates up to z ~ 120 
km, beyond which Pederson conductivity dominates Hall conductivity and penetration becomes 
mainly diffusive. 

For night time conditions (Fig. 0]), the currents and electromagnetic fields penetrate deeper up 
to z ~ 120 km by t = 0.1 sec. and vertical scale length of the structures increases. This is due 
to the decreased (by a factor of 5) conductivities which increase the wavelength and phase speed 
of the waves as A, oc 1/&h- The amplitudes of J x and J y decrease but not by factor of 5 as 
electric field amplitudes increase. In the relaxation phase, the structures have penetrated into the 
diffusive F-layer (above 120 km) up to the upper boundary of the simulation domain. The helicon 
waves with a typical wavelength of ~ 300 km propagate laterally in both directions. These waves 
are concentrated around the Hall conductivity maximum z = 110 km, which seems to guide the 
helicon waves. 

Fig. [5] shows results for continuous wave antenna field of frequency 10 Hz and 9 = 5° in the 
steady oscillatory state at t = 1.5 sec. This figure shows features similar to the case of the pulsed 
antenna field, viz., vertical wave penetration with decreasing wavelength and deeper penetration 
with larger wavelength for night time conductivities as compared to day time conductivities. How- 
ever, the vertical scale length is relatively smaller than that in the case of pulsed antenna field. 
The reason is that in the case of pulsed antenna field the Fourier spectrum of the antenna field 
consists of frequencies smaller than the frequency of the continuous wave antenna field. It can be 
seen from the wave dispersion relation that A oc 1/d which gives larger wavelengths corresponding 
to the lower frequencies in the pulsed antenna field case. 



A. ^-dependence 

Now we present results by varying 9 in order to see how the features observed for a single 
value of 9 = 5° change. The vertical profiles (along line x = 0) of the current density and electric 
field components are shown for various values of 9 in Figs. [6] and [7] (pulsed antenna field), and 
in Figs. [8] and [9] (continuous wave antenna field). The first thing that can immediately be noted 
from these figures is that the vertical scale length of the structures increases with 9; the structures 
become broader for larger values of 9. Consequently, the peak of the structures shift away from 
the boundary (z = 90 km) and deeper pentration takes place with increasing values of 9. The 
magnitudes of J y and J z decrease with 9 near the lower boundary. The vertical scale increases 
rapidly for up to 9 = 20°, but it is similar for 9 = 60° and 90° indicating asymptotic scale for large 
value of 9. For 9 = 60° and 90°, the profiles of J x and J y almost follow each other while profiles 
of J z have similar scale but differ in magnitude. The profiles of E x and E y also almost follow each 
other for 9 = 60° and 90°. The amplitudes of the spatial oscillations of E x and E y increase with 
9, while that of E z decreases very rapidly. For example, near the lower boundary z = 90 km, for 
pulsed antenna field and day time conditions E y increases from ~ —.05 mV/m (for 9 = 0°) to 
~ —0.1 mV/m (for 9 = 90°) and for continuous wave antenna field and day time conditions, it 
increases from ~ —.05 mV/m (for 9 = 0°) to ~ —.5 mV/m (for 9 = 90°). On the other hand for 
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FIG. 5. The magnetic field components B Xl By and B z (nT), the current density components j x -, jy and j z 
(nA/m 2 ) and the electric field components E x , E y and E z (mV/m) for continuous wave antenna field and 
9 = 5° in the steady oscillatory state (t = 1.5 sec). The left and right columns are for day and night time 
conditions respectively. 



pulsed antenna field and day time conditions E z decreases drastically from ~ 1.5 mV/m for 9 = 0° 
to very small values ~ 10 -4 mV/m for 9 = 90° and for continuous wave antenna field and day time 
conditions from ~ 3 mV/m for 9 = 0° to ~ 10 -4 mV/m for 9 = 90°. 

In order to understand the ^-dependence, we have derived a local dispersion relation from 
Eqs. (|13p - (|15p assuming uniform conductivities. Many of the observed features can be explained 
qualitatively from such a dispersion relation displayed below. 

p\u? + w/i [fci(po + Pp) + 2fcjjpp] = (k 2 ± + fcjj)[fcipppo + kf^pl + Ph)}, (19) 

where oj is frequency, k± = k z cos6 — k x sin9 and fen = k x cos9 + k z s\n9 (k x and k z are wave 
numbers along x and z directions respectively). In the collisionless limit (po — > 0,pp — > 0), this 
dispersion relation reduces to the whistler dispersion relation pquj = kk^p^- We have solved the 
dispersion relation (1191) numerically for k z for the given values of ui = 20tt sec. -1 , k x = 7r/200 
km -1 and values of po,pp and pn at z = 110 km where the Hall conductivity is maximum. The 
dispersion relation (|19j) is a fourth order polynomial in k z and hence its numerical solution gives 
four roots for k z . The two of the four roots are the complex conjugate of the other two. We 
discard the roots with negative imaginary parts as they give the solution exponentially growing in 
z. The real and imaginary parts of the other two roots as functions of 6 are plotted in top (root 
1) and bottom (root 2) panels of Fig. [lUJ As 6 increases, the real and imaginary parts of both the 
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FIG. 6. Vertical profiles of the current density components along the line cc = 0att = 0.1 sec. for pulsed 
antenna field and various values of 9. The left and right columns are for the day and night time conditions 
respectively. 



roots approach asymptotic values for large values of 9. For root 1, real part of k z drops down to 
very small values while imaginary part ~ 0.5 km -1 for large values of 9. This root has very long 
wavelength but decays in a small distance ~ 2 km. For root 2, imaginary part of k z attains very 
small values while real part ~ 0.4 km -1 for large values of 9, giving a large decaying distance and a 
wavelength ~ 15 km. This root will be visible in the simulations. The real part for root 2 decreases 
rapidly up to 9 ~ 20° and then slowly up to 9 = 90°. This is consistent with the simulations in 
which vertical scale length increases significantly up to 9 = 20° but there is not much difference in 
the vertical scales for large values of 9, e.g., scale lengths for 9 = 60° and 90° in Figs. [6]l9j 



B. A current distribution fit for J y 

It has been shown by Park [13] that for the purpose of the calculation of the electromagnetic 
fields produced by a current distribution J y ~ a/(x 2 + a 2 ) above the conducting Earth, J y can 
be replaced by a line current along y raised above its actual height by the half width (a) of the 
distribution. For this purpose, we integrate along z the two dimensional current distribution J y 
obtained from the simulations and fit a distribution ~ a/(x 2 + a 2 ) to the integrated current. The 
z-integrated current J y (x, t) = J J y (x, z, t)dz (mA/m) is shown in Fig. fTTh for 9 = 0°, day time and 
continuous wave antenna field. For a given time, J y has an x-profile which has its peak at x = 0. 
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FIG. 7. Vertical profiles of the electric field components along the line x = at t = 0.1 sec. for pulsed 
antenna field and various values of 9. The left and right columns are for the day and night time conditions 
respectively. 



This peak value J y {x = 0, t) oscillates in time with amplitude 1.4 mA/m in steady oscillatory 
state, shown in Fig lllb which also shows the fit p(t) for J y (x = 0,i). This fit p{t) has the same 
functional form as the continuous wave antenna field except the amplitude in the steady oscillatory 
state, which is 1.4 mA/m for p(t). Using p(t), J y (x,t) is fitted with the function, 

J^(x,t) = 4^%, (20) 

where a = 100 km. The plots of J y (x,t) and J y lt (x,t) at two different times (when p(t) is at its 
positive and negative peaks ) are shown in Fig. [Tib . The total current I y flowing along y obtained 
by integrating J y (x,t) and J y lt {x,t) along x, shown in Fig. [TTb l. oscillates with amplitudes 310 
Amps and 440 Amps respectively. The mismatch between the total currents is due to the slight 
mismatch between J y (x,t) and J y lt (x,t) for x > 100 km. 

The z-integrated current shown in Fig. [11] is for 8 = 0°, day time conditions and continuous 
wave antenna field. However the peak value of p(t) (1.4 mA/m), the fit J y lt and the total current 
I y do not change at all (except the time dependence of p(t) which is different for pulsed and 
continuous wave antenna field) on changing the simulation parameters, viz., day and night time 
conductivities, pulsed and continuous wave antenna field, and values of 6. This can be understood 



13 



ecr 



-500 

90 



Day Time 











J (riA/m 2 ) 



200 
1 OO 
O 

-1 OO 



90 u 

Night Time 







^-A J (nA/m 2 ) 



1 OO 



1 20 



1 40 



90 



1 1 O 



1 30 



1 50 




1 40 



200 
1 OO 
O 

-1 OO 

-200 
9 











J (nA/m 2 ) 



1 1 o 



1 30 



1 50 




1 40 









Vy \ / J z (nA/m 2 ) 



90 



1 1 O 



(km) 



1 30 



1 50 



FIG. 8. Vertical profiles of the current density components along the line x = 0att=1.5 sec. for 
continuous wave antenna field and various values of 9. The left and right columns are for the day and night 
time conditions respectively. 



as follows. Since d/dx « d/dz and B z is small as compared to B x , J y can be approximated as, 

i dB x 

fio dz 

On integrating from lower boundary at z = 90 km to the upper boundary we get, 

Jy(x,t) = I J y (x,z,t)dz = — B x (x, z = 90km, t), 

using B x = at the upper boundary. Due to the magnetic shielding, the value of B x ~ 1.8 nT 
at x = and z = 90 km (when antenna field is at its maximum) irrespective of the simulation 
parameters. This gives the peak value of J y {x = 0,t) ~ 1.4 mA/m, which is same as in Fig. Illb . 
The fit of Jy along x is basically related to the x-dependence of the x-component of the antenna 
field (B Xtan t) for a given z (see Eq. [8|). Since the peak value 1.4 mA/m and x-dependence of the 
B x ,ant remain same for all the simulation parameters, the total current I y also remains same. 



IV. SUMMARY 

We have investigated the interaction of the ELF electromagnetic fields generated by an antenna 
placed at the Earth's surface with the E-region of the Earth's ionosphere for various parameters, 
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FIG. 9. Vertical profiles of the electric field components along the line x = at t = 1.5 sec. for continuous 
wave antenna field and various values of 9. The left and right columns are for the day and night time 
conditions respectively. 



viz., day and night time conductivities, pulsed and continuous wave antenna field, and different 
values of 9. The interaction leads to the generation of the horizontal currents J x and J y with 
magnitudes up to few 100 nA/m 2 and vertical currents J z with magnitude up to 7 nA/m 2 in 
the E-region depending upon various simulation parameters. The associated electric fields have 
magnitude up to E x ~ 0.6 mV/m, E y ~ 1 mV/m and E z ~ 8 mV/m. The wave penetration 
with a typical wavelength of the order of ~ 10 km of the currents and fields in to the deeper 
ionospheric layers up to z ~ 120 km takes place due to the dominance of the Hall conductivity over 
the Pederson conductivity in the region 90 km < z < 120 km and penetration becomes diffusive for 
z > 120 km. The scale length A of the vertical penetration decreases with z up to z = 110 km as 
A oc 1/oh and an increases up to z = 110 km. The increase in the wave speed due to the reduced 
conductivities during night time leads to the deeper penetration. For continuous wave antenna 
field, the vertical scale length is smaller than that in the case of pulsed antenna field because of 
the presence of lower frequencies in the later case. The vertical scale length increases rapidly with 
9 up to 9 = 20° and then slowly to approach an asymptotic value for large values of 6. This leads 
to the deeper penetration with increasing value of 9. The magnitudes of J y and J z near the lower 
boundary decrease with 9. The vertical electric field E z decreases drastically with 9, e.g., from 1.5 
mV/m [9 = 0°) to 10~ 4 mV/m (9 = 90°) for pulsed antenna field and day time conditions. The 
z-integrated current along y is fitted with a current distribution a 2 p(t)/(x 2 + a 2 ) where a = 100 
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FIG. 10. ^-dependence of the real and imaginary parts of the vertical wave number k z obtained from the 
solution of the dispersion relation (p~9|) for given values of cj = 2Qtt s -1 , k x — 7r/200 km" 1 and values of 
pp, pn and pq at z — 110 km. 



km and p(t) is the peak value of the distribution. Such a current distribution can be replaced by 
a line current raised above its actual height by the half width (a) of the current distribution for 
the purpose of the calculation of radiation 171 ] . The maximum total current along y (310 Amps) 
remains same for various simulation parameters due to the magnetic shielding which gives B x ~ 1.8 
nT at x = and z = 90 km for all the simulation parameters. 
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